(De) confinement of supercurrent in Z2 topological Insulators 
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Abstract 

It is shown that a Z2 topological system with C/cm(l) x Uz{l) (electromagnetic and spin) gauge 
symmetries shows superconductivity and quantized spin Hall effect simultaneously. When Uz{^) is 
broken, a dissipationless electric current is still possible to flow locally but net charge transfer is 
absent, i.e., current is confined. In the Kane-Mele model for graphene, this confining-deconfining 
transition of the supercurrent is driven by the Rashba spin-orbit interaction, which breaks Uz{l)- 

PACS numbers: 



* Present Address: Institute for Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro, Tokyo 
153-8505 



1 



I. INTRODUCTION 



Dissipationless transport in a certain class of band insulators is often called "topological" . 
We..-U„ow„ examples are integer and tacfona. quantum Hal. effectsQ. specific type of 
anomalous Hall effect [2|, etc. Dissipationless electric current in these examples is associated 
with time-reversal symmetry (TRS) breaking. However, topological transport for spin is not 
accompanied by breaking of TRS. The Kane-Mele (KM) model for graphenejs, 4| preserves 
TRS, and realizes a quantized spin Hall (QSH) effect fl. Unfortunately, quantization of spin 
Hall conductivity (SHC) is not robust against perturbations which break the residual spin- 
rotation symmetry f/^(l), e.g., Rashba spin-orbit interaction (SOI) 6|. In the presence of 
such perturbations, the role of spin Chern number, characterizing the QSH state is replaced 



by a Z2 topological invariant 



aa. 



Physically, the same role is played by a pair of helical 



ed ge m odes [j, |8|. Z2 classification of topological insulators is extended to three dimensions 

0,13 • 



Z2 topological insulators show unique electromagnetic properties ll|. In contrast to Ref. 
which takes into account only topologically robust quantities protected by Z2 invari- 
ants, this paper reveals the role of non-robust QSH effect in electromagnetic responses of 
Z2 topological insulator. We show that QSH effect results in/from massive electrodynamics 
characteristic to superconductivity. This mechanism of superconductivity is similar to the 
arguments in a TRS system 12| with "topological order" characterized by the ground state 
degeneracy 13|] • On a torus geometry, the superconducting phase corresponds to deconfine- 
mentoi a dissipationless current (non-contractible current loops exist; see Fig. [T]). Whereas 
a more generic Z2 phase {\r ^ case in KM model) corresponds to "confining" state, where 
the current can flow locally (contractible current loops exist; see. Fig. [1]), but no net cur- 
rent flows around the torus. Thus, we find that is a "switch" for dissipationless charge 
transfer. 

Let us consider the following Hamiltonian of the KM model for graphene 0, 4| in the 
continuum limit, l-i = 'Ho + "Ha + 'Hr, which consists of the following three elements: 



Ho 
-Ha 



X 



R . 



(1) 
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where and Dy are covariant derivatives defined below in Eq. (|3]). Physical origins of each 
element composing "H are the following: I-Lq describes idealized graphene in the massless 
limit. gives a topological mass gap A, which stems from the intrinsic part of SOI. I-Lr is 
the Rashba term, induced by an electric field applied perpendicularly to the graphene layer, 
and breaks the spin-rotational symmetry, i.e., the conservation of Sz- Note that there are 
three species of (pseudo) spins in the model: A-B sublattice spin a, valley (K-K') spin r, 
and real spin s, where, e.g., a = {crx,cry,crz) represent Pauli matrices. 
Our total Lagrangian density reads, 

£ = ^t(^zDo-^)^ + ^(^'-52), (2) 

in which we take into account the dynamics of electromagnetic fields (Maxwell term). 
represents the electron field. The electromagnetic field (/i = 0, x, y) couples to ^ 
through a covariant derivative in C defined by 

iD^ = idf, - eA^ + ya^. (3) 

Here we have introduced a gauge field a^, which couples to the spin degrees of freedom of \E'. 
The spin gauge field a^, thus introduced, plays the central role in the following discussion. 
Let us consider a functional integral with an additional term a^J^ in Eq. ([2]), i.e.. 

Since the (2+l)-dimensional ((2+l)D) spin current is defined by = —dC/da^, the La- 
grangian is rewritten as £ = Ca^=o — a^J^p. Then, integrating out in Eq. (j4]), one 
finds 

e'^''^^^'''^'^^ = j V^>^V^>5{J% - J^p)e^^'^'^'^"''-°, (5) 

i.e., a constraint, J^'^ = J^, is imposed on the spin current. Hence, Scs[A^, is regarded 
as an effective action for the electromagnetic fields under a given (2+l)D spin current Jj^^. 
Eq. (|5]) thus clarifies the physical meaning of introducing the gauge field a^. Using gauge 
field is a standard prescription in field theory for implementing some constraints on physical 



quantities 13| . We will discuss elsewhere that we can also introduce Cq and its integration to 
describe the electron correlation via Stratonovich-Hubbard transformation jl^ . The result 
in the strong correlation limit is, in essence, equivalent to one shown in this paper. 
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In the following, we derive the explicit form of Scs[A^, Jtp]y discuss unique electromag- 
netic response properties of the QSH insulator. We consider, separately, (i) s^-conserving 
{Xr = 0), and (ii) s^-non-conserving (A/j 7^ 0) cases. 



II. s^-CONSERVING CASE: "SUPERCONDUCTIVITY" AND QSH EFFECT 



When Xr = 0, the system has a local f/cm(l) x Uz{l) symmetry, i.e., the Lagrangian 
density is invariant, independently, under \1/ — )■ e~*^^\l/,/l^ — )■ 
^ _^ Qis,'fi/2^^ _s. + g^^p_ Integrating out ^ in Eq. (jl]). 



dfj_6, and under 



one finds, in the Gaussian approximation 



-p^ ^iS,BlA'^,a'^,Jtp] 



(6) 



2 2/i 



We have introduced the field strength, Cj = SoOj — ditto, and 6 = eijdiaj = x,y), 
associated with the gauge field a^. Note that the dielectric constant e and the magnetic 
permeability fi are renormalized: e = 1 + l/GvrA, /i = (1 + f|,/67rA)~^. Here, we use the 
unit in which vacuum values of these constants are unity: eo = /io = 1- 

The first terms (a A dA + A A da) on the r.h.s. of Eq. (171) are the so-called BP term 

n tt 

|15l |. Like the Chern-Simons (CS) term {A A dA, or a A (ia)[16|, these terms are topological, 
i.e. independent of the metric of a manifold. The CS term breaks TRS and is not induced 
in Eq. ([7]). The BP term preserves TRS. Recall that TRS operation transforms Aq — )■ Aq, 
A — —A, whereas, ao — )■ — oo? « — ?■ a. 

The coefficient in front of the BP term in Eq. ([7]) is precisely the SHC a^y (up to a 
symmetrization factor 1/2); the latter obtained by Kubo formula. The physical impulication 
of BP term is, thus, the spin Hall effect. Note that SHC per a valley is (e/47r)sgn(A), which is 
related to the winding number of Permion propagator in (2+l)D momentum representation 



17|. If one considers a Bloch system, i.e., a momentum space compactified on a two- 



torus, one finds, after frequency integration [l8|, a^y = ^Csc, where Csc is the spin Chern 



integer [6| (See, also Refs. 19|, |20|). This quantization is assured whenever Sz is conserved. 



We can see the bulk/edge correspondence as a consequence of the presence of t/em(l) x 
Uz{l) gauge symmetry and quantum anomaly, which is an analogy of the field-theoretic 



(7) 



4 



argument of the edge state in quantum Hall system 2l| (see, also 22|). The number of 
helical edge modes is counted by Csc, equivalently by SHC cr*^. This fact has been already 
verified by Ref. 6| numerically. 

Note also that a kinetic term of (Maxwell term associated with) is also induced in 
Eq. d?]). Its coefficient in front is e — 1 (in contrast to e for A^), indicating that the term 
arises due to the renormalization of dielectric constant. The gauge field thus acquires its 



own dynamics. Some related arguments have been proposed 



12 



23| 



15|, analogous to the CS termjl6l.l26|. 



The BF term gives the massive electrodynamics 
We proceed to integrate out in Eq. (E]) using the Lorentz gauge. Then, we see that, for 
0, 



(TX 



(A--.^A-) 



52 
2/i 



(8) 



The superscript T denotes the transverse component: Al^ = ~ 
(9^/n) [(l/vpydoAo — djAj], where we introduced a d'Alembert operator (d'Alembertian) 
suitable for graphene: □ = (l/t'i?)^9o — df. The presence of Aq"^ and Af^ terms in Eq. 
([8]) dictates that the electromagnetic field becomes indeed massive. These terms are 
known to be induced in the effective action for superconductivity, after integrating out the 
Nambu-Goldstone mode 2J1 (see also Ref. 2511'). 

Note that (possible) superconductivity in topological insulators mentioned in Ref. [27[ is 
different from ours. The former is predicted outside the QSH phase. In our scenario, super- 
conductivity arises in the QSH phase, and is shown to be indispensable for the occurrence 
of QSH effect. 

Let us examine the QSH effect. In the presence of finite spin current, J^^ 7^ 0, an 

C\n — — _ 

additional term 



, ^ffA^^i^Jg'; + C(Jsf ), is in the effective action (ED. Suppose that we 



apply a static spin current in the y-direction: J^^ 
equation of motion for E reads. 



J: 



sp 



dx 



:eEJx) 



^-^E.{x) + ^^Jl{x). 



0. Then, the 



(9) 



In the long- wavelength limit, the second order derivative term on the l.h.s. of Eq. ^ is 
irrelevant, leading to the QSH effect. The presence of o"^^ term coming from the gauge mass 



term is essential. The equation o 



for a thin-film superconductor 
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motion for the magnetic field B is analogous to the one 



III. s^-NON-CONSERVING CASE: CONFINEMENT OF DISSIPATIONLESS 
CURRENT 

When A/j 7^ 0, the SHC is not quantized, and loses its topological meaning {g]. However, it 
still plays an important role for the electromagnetic response. The action S'cff[^^, a^-, in 

2 

Eq. ([7]) obtained after fermion integration is modified and a mass term of a^, ^(oq — Vpof), 
appears as a direct consequence of the breakdown of local Uz{l) symmetry. Here, = 
(A|j/87r|A|) + 0(A|j,). SHC, and other coefficients e and fi in Eq. ([7]) receive corrections in 
the order of C(A^). Thus, SHC is no longer quantized. By integrating out for = 0, 
we obtain 



f AT ^ AT AT --^^ 



-A 



Xe-l)vl\ °n + m2 ° n + ml" J ' 2 2/i 



(10) 



If we introduce Jgp 7^ 0, the spin Hall effect occurs with non-quantized SHC. 

We first introduce the notion of confinement of dissipationless current. We consider a 
system on a two-torus to avoid unnecessary complications due to the presence of an edge. 
The Fourier transform of static charge current density in Eq. (iTOll is 



J. = - '\ ,1 , A,, (11) 

which is dissipationless, since it is driven by a static vector potential that never yields a 
voltage drop. Eq. f|TT]) also indicates that current density j{x) can be finite, on the other 
hand, the net current at an arbitrary cross section of the torus is zero, i.e., current is 



confined 



311 . The latter can be shown by using the current conservation relation and the 



fact that jq=o = 0. We note that jq=o is equal to the integral of current density over the 
entire surface of the torus. Absence of net current is a usual insulating behavior. Unusual 
fact in this confining state is that a finite current can fiow locally. Such a local current 
circulates around a contractible closed contour (see Fig. [I](a)). 



IV. DECONFINING (SUPERCONDUCTING) PHASE REVISITED 

We have seen that spin-rotation symmetry breaking yields the confinement of dissipa- 
tionless charge current. On the other hand, in the spin-rotation symmetric limit, — 0, 
Eq. ffTTl) reduces to, jg = —-^Aq, i.e., identical to the supercurrent obtained from Eq. (|8]). 
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(a) Confining state (b) Deconfining state 
(superconducting) 




FIG. 1: Schematics of possible current loops in (a) confining state, and (b) deconfining state for 
dissipationless electric current. 



TABLE I: Topology of possible current loops on the torus. 





current flow along a 




contractible loop 


non-contractible loop 


QSH phase 


NO 


YES 


(Aif^ = 0) 




— )• deconfining phase 


Z2 phase 


YES 


NO 


(Ai? + 0) 




—7- confining phase 



This means that j,j=o can be finite, allowing for a finite net current circulating around the 
torus, i.e., "deconfinement" is accomplished. 

In the deconfining state, non-contractible current loops become possible (see Fig. [I](b)). 
On the other hand, contractible loops are forbidden because of the superconducting nature of 
the state (perfect diamagnetism) . This contrasting behavior is summarized in Table HI The 
topology of possible current loops thus plays the role of "order parameter" which classifies 
the confining and deconfining phases. This fact reminds us the role of Wilson loop in the 



lattice gauge field theory 



29|. 



This confinement- deconfinement transition is expected in graphene as the Rashba cou- 
pling Aj^ is switched on by an electric field perpendicular to the system. It would be interest- 
ing if the transition can be applied to the development of switching devices for dissipationless 
charge transfer. It should be noted that confinement occurs even for the infinitesimally small 
spin- rotation symmetry breaking (see Eq. f|TT]) ). Then, a perfect-spin-conserving system 
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would be needed to observe the deconfining (superconducting) state. 



V. HGTE/CDTE QUANTUM WELL 

Although we have used the KM model for demonstration, the basic ideas presented in 
this paper will be applicable to other QSH insulators. The present result with spin-rotation 
symmetry would be relevant to HgTe/CdTe quantum well js], [sO], which has a large band 
gap (~100K) in bulk spectrum assisting the experimental observation. In this system, we 
can define a conserved pseudo spin associated with the band index, although real spin is 
not conserved. In parallel, we can see that the superconductivity is induced by integrating 
out Fermion and pseudo spin gauge field. Although protected from the thermal activations 
by the large band gap, this superconductivity is fragile against the pseudo spin-rotation 
symmetry breaking perturbation. 



VI. SUMMARY 



It is shown that the KM model in s^-conserving case shows superconductivity and quan- 
tized spin Hall effect simultaneously. When s^-conservation is broken by the Rashba spin- 
orbit interaction, a dissipationless electric current is still possible to flow locally but net 
charge transfer is absent, i.e., current is confined. 

We will discuss elsewhere that we can also introduce the 0-th component of the spin 
gauge fleld and its integration to describe the electron correlation via Stratonovich- Hubbard 
transformation. The result in the strong correlation limit is essentially equivalent to one 
which has been shown in this paper. 
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